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ABSTRACT 


This  investigation  attempts  to  define  the  flexural  and 
lateral-torsional  buckling  strengths  of  double  angle  struts. 

CSA-G40.12  steel  was  used  to  establish  the  basic  material 
properties.  Residual  stress  measurements  were  performed  on  specimens 
cut  from  five  different  angle  shapes.  The  results  are  reported 
herein  along  with  the  results  of  tension  tests  on  comparison  specimens. 

The  residual  stress  distribution  was  idealized  and  used  as 
input  for  a  computer  program  which  computed  the  flexural  and  lateral - 
torsional  buckling  strengths  based  on  the  tangent  modulus  concept. 
Column  curves  were  developed  for  several  different  column  sections 
to  investigate  the  influence  of  member  asymmetry,  leg  thickness, 
angle  separation,  yield  stress  level  and  residual  stress  distribu- 
ti  on . 

For  the  double  angle  members,  the  lateral -torsional  buckling 
strength  coincided  (within  2%)  with  the  (fictitious)  y-y  axis  buck¬ 
ling  strength  for  all  cases  investigated.  Thus  the  traditional  de¬ 
sign  which  ignores  the  possibility  of  lateral-torsional  buckling  may 
be  used  with  confidence  for  these  members. 
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CHAPTER  I 

INTRODUCTION 

Double  angles  are  commonly  used  in  building  structures  as 

compression  members  in  trusses  and  as  bracing  members  to  ensure  the 

stability  of  the  structure.  These  members  are  normally  considered  to 

be  loaded  axially,  as  the  secondary  bending  moments,  induced  by  the 

1 5 

connections,  are  small 

A  typical  building  panel  is  shown  in  Figure  1.1.  The  double 

angles  are  used  as  diagonal  bracing  members  and  may  be  designed  to  act 

only  in  tension,  or  to  resist  both  tensile  and  compressive  forces.  The 

angles  forming  a  member  are  separated  by  the  thickness  of  the  gusset 

plate  (to  facilitate  connection)  and  are  either  bolted  or  welded  to 

the  connection  plate.  A  typical  bolted  connection  is  shown  in  Figure 

1.1(b)  and  the  corresponding  welded  detail  in  Figure  1.1(c).  The  angles 

are  connected  at  intervals  along  their  lengths  so  that  the  two  angles 

act  as  a  single  unit.  The  end  restraint  offered  by  the  usual  type  of 
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connection  is  small  .  Thus  in  design,  it  is  assumed  that  the  member 
acts  as  a  pin-ended  column  having  an  effective  length  equal  to  the 
distance  between  the  extremities  of  the  two  connections. 

The  range  of  angles  sizes  used  in  building  structures  generally 
varies  from  the  very  small  (3  x  2y  x  yg-")  to  the  larger  size  (8  x  8  x  lg"), 
with  a  selection  of  equal  and  unequal  leg  sizes  available  between  these 
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limits.  The  length  of  the  angle  leg  governs  the  location  of  the  bolts 

in  the  connection.  Standard  gauge  distances  are  used  as  shown  in 

Table  1.1,  with  g  representing  the  gauge  distance  for  a  single  line 

of  bolts  and  g-j  and  g^  denoting  the  gauge  distance  if  two  lines  of 

bolts  are  used.  The  use  of  welded  connections  as  shown  in  Figure 

1.1(c)  is  also  common  with  the  leg  thickness  governing  the  size  and 

5 

strength  of  the  weld  . 

The  ultimate  strength  of  a  compression  member  is  usually 
assumed  to  be  equal  to  its  buckling  strength;  defined  as  the  load  at 
which  the  member  can  exist  in  two  adjacent  equilibrium  positions. 

Figure  1.2  shows  a  schematic  column  curve  which  plots  the  applied 
stress  corresponding  to  the  buckling  load,  a,  non-dimensional ized  by 
the  yield  stress,  o  ,  versus  the  slenderness  ratio  (L/r) ,  where  L  re- 
presents  the  member  length  and  r  the  radius  of  gyration. 

The  branch  E-D-B-C  of  the  curve  represents  the  relationship 
between  the  critical  applied  stress,  a,  and  the  slenderness  ratio  when 
buckling  occurs  before  yielding  of  any  fibre  of  the  cross-section. 

For  an  el astic-perfectly  plastic  material  this  curve  would  remain  valid 
until  the  value  of  o  reaches  the  yield  stress,  a  ,  of  the  material. 

sJ 

For  smaller  values  of  L/r,  the  critical  stress  would  remain  at  a  thus 
defining  the  branch  A-D-B-C. 

Residual  stresses  caused  by  the  hot  rolling  and  cooling  process, 
cause  sections  of  the  cross-section  to  yield  at  stresses  below  o  .  This 
produces  a  deviation  from  the  el asti c-perfectly  plastic  curve,  at  a 
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stress  level  equal  to  yield  stress  minus  the  maximum  residual  stress, 
a  ,  as  shown  by  the  curve  A-B-C.  As  the  load  is  increased  parts  of  the 
cross-section  yield  before  buckling  occurs.  For  smaller  values  of  L/r, 
larger  portions  of  the  cross-section  yield  before  buckling  until  finally 
the  entire  cross-section  is  yielded  at  a/a  =  1.0. 

<y 

Buckling  of  doubly  symmetric  sections  may  occur  in  any  one 
of  three  independent  modes;  flexural  buckling  about  either  of  the  two 
principal  axes  or  torsional  buckling.  Double  angle  struts  have  only  a 
single  axis  of  symmetry,  however,  and  thus  there  exist  only  two  in¬ 
dependent  modes;  flexural  buckling  about  a  principal  axis,  or  a  com¬ 
bination  of  torsional  buckling  and  flexural  buckling  about  the  other 
principal  axis  (lateral-torsional  buckling)  as  shown  in  Figure  1.4. 

The  present  investigation  attempts  to  define  both  the  flexural 
and  lateral-torsional  buckling  strengths  of  double  angle  struts.  The 
results  will  be  compared  with  design  provisions  with  the  aim  of  deter¬ 
mining  the  factors  of  safety  inherent  in  the  various  procedures. 
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Figure  1.2  Buckling  Strength  of  Compression  Member 
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Figure  1.3 


Modes  of  Buckling  for  Double  Angle  Struts 
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Table  1.1  Standard  Gauge  Distance  for  Angles 


CHAPTER  II 


PREVIOUS  INVESTIGATIONS 

The  solution  to  the  differential  equation,  which  expresses 
the  equilibrium  of  the  member  in  the  deformed  position  (Figure  2.1) 
was  first  obtained  by  Euler  in  1759^,  The  solution  defines  the  load 
at  which  a  slender,  axially  loaded  column  of  uniform  cross-section  can 
exist  in  either  a  straight  or  slightly  deformed  configuration.  The 
solution  is  limited  to  slender  columns  in  which  the  strain  at  the 
instant  of  buckling  does  not  exceed  the  yield  strain. 

Most  practical  columns  do  not  fall  within  the  range  of 
applicability  of  the  Euler  solution  and  in  1889  Engesser^  developed 
the  first  inelastic  budding  theory.  This,  the  original  tangent  modul 
theory,  assumed  that  the  column  remained  straight  until  the  buckling 
load  was  attained  and  that  the  modulus  of  elasticity  associated  with 
the  buckling  motion  was  the  tangent  modulus  at  the  critical  stress 
level . 

In  1895  Considere  and  Jasinski  challenged  this  theory  on  the 
basis  that  lateral  deflection  would  produce  additional  loading  (due 
to  bending)  on  the  concave  side  and  a  corresponding  unloading  on  the 
convex  side.  In  the  light  of  this  criticism  Engesser  modified  his 
original  thinking  and  postulated  the  reduced  modulus  theory,  wherein 
the  tangent  modulus  was  applied  to  only  the  loading  portion  of  the 
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cross-section  and  the  original  elastic  modulus  was  used  to  determine  the 
stiffness  of  the  unloading  portion.  Tests  conducted  by  von  Karman  in 
1910,  however,  agreed  well  with  the  original  tangent  modulus  theory^. 

The  difference  between  the  two  concepts  remained  unresolved 
until  1947.  From  the  results  of  careful  model  tests  Shanley  showed 
that  an  initially  straight  column  will  buckle  at  the  tangent  modulus 
load,  then  continue  to  deform  with  increasing  load.  Shanley  concluded 
that  for  an  initially  straight  column,  the  tangent  modulus  load  is  a 
lower  bound  on  the  ultimate  strength  of  the  column  and  that  the  re¬ 
duced  modulus  load  gives  an  upper  bound  for  the  ultimate  strength,  if 
the  column  is  temporarily  supported  beyond  the  tangent  modulus  load. 

The  ultimate  strength  of  an  actual  column  lies  between  the  tangent 
modulus  and  reduced  modulus  load,  with  the  tangent  modulus  being  a 
reasonable  estimate  of  the  strength  for  structural  steel  members. 

The  results  of  Shanley 's  work  led  to  a  reinstatement  of  the  tangent 
modulus  theory,  and  to  additional  applications  of  the  basic  concept,  in 
that  strain  reversals  of  previously  yielded  fibres  are  ignored  in 
computing  the  buckling  strength. 

Residual  stresses  affect  column  strength  in  that  the  "ef¬ 
fective  "stress-strain  curve  for  the  complete  cross-section  is  en¬ 
tirely  different  from  the  stress-strain  curve  obtained  from  a  coupon 
taken  from  the  material.  In  preparing  tensile  coupons  for  testing, 
residual  stresses  are  at  least  partially  released  and  the  stress- 
strain  curve  approaches  the  idealized  elastic-plastic  curve.  To  assess 
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the  magnitude  of  the  residual  stresses,  stub  column  tests  on  the  com- 

g 

plete  cross-section  have  been  used  .  In  a  stub  column  the  residual 
stresses  cause  individual  fibres  to  yield  well  below  the  attainment 
of  the  yield  load  for  the  specimen.  The  yielded  fibres  offer  no  re¬ 
sistance  to  additional  loading  and  thus  the  load-deflection  (stress- 
strain)  curve  deviates  from  the  original  (elastic)  straight  line. 

The  stress  at  the  point  of  deviation  then  can  be  used  as  a  crude 
indication  of  the  maximum  compressive  residual  stress.  The  alter¬ 
native  to  this  procedure  is  to  determine  the  distribution  of  residual 
stresses  using  a  sectioning  procedure.  The  use  of  the  tangent  modulus 

theory  is  equivalent  to  computing  the  flexural  stiffness  of  the  member 

4 

based  on  the  elastic  core  of  the  cross-section  at  buckling  .  The 
critical  stress  level  can  then  be  formulated  as 


(2.1) 


where  E,  represents  the  modulus  of  elasticity  and  I  is  the  moment  of 

inertia  of  the  elastic  core.  The  cross-sectional  area  is  represented 

by  A  while  L  denotes  the  length  of  the  column. 

The  basic  equations  for  both  flexural  and  torsional  buckling 

have  been  derived  for  doubly  symmetric  sections  and  have  been  solved 

for  the  buckling  strength  of  the  member  in  both  the  elastic  and  inelastic 

1 

ranges.  Lee  et  a]_  have  reported  the  results  of  an  analytical  study 
of  the  buckling  strength  of  wide  flange  sections.  In  this  study  the 
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influence  of  the  residual  stress  distribution  on  the  torsional  buckling 
strength  is  reported.  Three  conditions  must  be  satisfied  by  the  re¬ 
sidual  stress  distribution  in  order  that  the  member  is  in  equilibrium 
before  the  load  is  applied.  The  net  axial  force  acting  on  the  member 
must  be  equal  to  zero;  this  implies  that: 


afdA  =  0  (2.2) 

A 

In  addition  the  bending  moments  about  the  x  and  y  axes  must  vanish; 
this  implies  that  both: 


and 


arxdA  =  0  (2.3) 

c 

A 


<* 

c^dA  =  0  (2.4) 

# 

A 


Lee  also  requires  that  the  residual  stress  distribution 
produce  no  twisting  moment  with  respect  to  the  shear  centre  for  any  small 
deformation.  Thus 


ar  <x02  +  OdA  =  0  (2-5) 

A 

where  xQ  and  yQ  represent  the  location  of  the  shear  centre.  If  this 
fourth  condition  is  satisfied,  there  will  be  no  reduction  in  the 
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torsional  buckling  load  due  to  residual  stresses,  except  that  reduction 

associated  with  the  gradual  deterioration  of  the  cross-section. 

1 2 

It  has  been  shown  that  Equation  2.5  is  not  a  necessary 
condition  for  equilibrium  of  the  residual  stresses  and  in  fact,  measured 
distributions  do  not  satisfy  this  condition.  This  implies  that  the 
presence  of  residual  stresses  will  reduce  the  torsional  buckling 
strength  somewhat  and  must  be  considered. 

The  deformations  at  the  instant  of  (flexural-torsional) 
buckling,  for  singly  symmetric  sections,  involve  a  translations, 
in  the  x-di recti  on  as  shown  in  Figure  2.2  and  a  rotati on, 0, about  the 
shear  centre.  Pure  flexural  buckling  of  the  same  section  involves 

g 

only  a  translation  in  the  y-di rection  v.  Chajes  and  Winter  have  sum¬ 
marized  the  three  equilibrium  equations  for  a  singly  symmetric  section 
in  the  deformed  condition  as: 


Ely  uIV  +  P(uH  +  y04>n)  =  0 

(2.6) 

El  vIV  +  P(vH)  =  0 

A 

(2.7) 

EIW  <t>1V  -  (GKt  +  KH11  +  PyQ  u11  =  0 

(2.8) 

where  I  and  I  are  the  moments  of  inertia  about  the  principal  axes, 
x  y 

I  represents  the  torsion  constant  for  warping  and  Kj  the  St.  Venant 
torsion  constant.  The  deformations  are  given  by  u,  v  and  <|>  (Figure  2.2) 
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and  the  primes  represent  differentiation  with  respect  to  the  co-ordinate 

measured  along  the  longitudinal  axis.  The  axial  load  P  is  given  by 

/  adA  and  K  is  defined  as  /  aa  dA  where  a  represents  the  distance  from 
A  A 

a  point  on  the  cross-section  to  the  shear  centre. 

In  the  derivations  of  equations  2.6  to  2.8  it  was  assumed 
that  the  deformations  were  small  and  that  the  material  was  elastic. 
Solutions  to  these  equations  lead  to  the  two  buckling  loads  for  the 
member.  For  a  column  with  ideally  pinned  boundary  conditions  at  both 
2=0  and  2  =  L; 


u  =  v  =  4)  =  0 


and 


u 1 '  =  v '  '  =4)"  =  0 


The  equations  result  in  the  following  solution: 


-  f  (p-p  )(p-px)(p-p*)  -  pV^-v  =  0 


where 


tt2EI, 

p  =  - • 

y  i2 


2CT 
TT  El 

P  = 
x  L2 


X 


it2EI 


P<t>  = 


-  =  [GKj  +  — 

K  L 


(2.9) 


(2.10) 

(2.11) 


and 


(2.12) 
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Equation  2.9  has  two  solutions,  P  =  P  which  represents  pure  flexural 

A 

buckling  and  P^  equal  to  the  lower  root  of  the  quadratic  resulting 
from  the  remainder  of  Equation  2.9.  This  load  represents  the  coupled 
flexural-torsional  buckling  load  and  is  always  below  either  P,  or  P 

(p  y 

The  lower  of  the  two  solutions  of  Equation  2.9  represents  the  critical 
buckling  load  for  the  member. 

The  solution  of  Equation  2.9  in  the  elastic  range  is  straight 
forward;  the  solution  of  Equation  2.9  in  the  inelastic  range  is  valid 
if  the  section  properties  used  in  Equations  2.10  to  2.12  are  based  on 
the  elastic  core  remaining  at  each  load  interval.  The  St.  Venant 

o 

torsion  resistance  is  based  on  the  original  section  .  The  extent 
of  the  remaining  elastic  core  remaining  at  each  load  level  depends  on 
the  residual  stress  distribution.  Residual  stresses  are  induced  in  a 
wide  flange  section  by  differential  cooling  and  plastic  flow  after  hot 
rolling.  The  flange  tips  cool  and  stiffen  first,  then  the  web  cools 
and  must  contract  against  the  previously  solidified  flanges.  This 
cooling  sequence  results  in  tensile  residual  stress  in  the  web  and 
compression  stresses  at  the  flange  tips.  Similar  residual  stress 
patterns  are  formed  in  angle  sections  as  a  result  of  differential 
cool ing. 

Maximum  values  of  residual  stresses  in  rolled  wide  flange 

14 

sections  have  been  previously  reported.  Beedle  and  Tall  report  the 
average  value  of  the  peak  compressive  residual  stress  as  0.3  of  the 
yield  stress  for  ASTM-A7  steel.  Peak  values  in  angle  sections  have  not 
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9 

been  investigated  to  any  extent.  However,  Galambos  reports  the  results 
of  two  stub  column  tests  on  single  angles,  giving  an  approximate  peak 
compressive  residual  stress  of  0.30  and  0.22  times  the  static  yield 

3 

stress.  Beedle  and  Tall  report  a  residual  stress  pattern  for  angles 

in  riveted  columns  (Figure  2.3)  with  a  peak  value  of  approximately  0.3 

1 8 

times  the  yield  stress.  O'Connor  measured  residual  stresses  in 
random  samples  of  angles  (Figure  2.4)  in  an  attempt  to  establish  the 
residual  stress  distribution.  The  stresses  reported  in  these  tests 
did  not  conform  to  any  consistent  pattern. 

Test  data  relating  the  strength  of  double  angle  struts  with 
or  without  residual  stress  considerations  are  almost  non-existent. 
Published  work  on  the  buckling  strength  of  other  singly  symmetric 
sections  in  the  inelastic  range  is  also  meagre. 

g 

Galambos  reported  the  results  of  tests  on  single  angles 
loaded  axially;  the  sections  deflected  in  the  principal  directions 
with  only  minor  twisting  deformations  observed.  Lenzen^  has  tested 
open  web  steel  joist  sections  with  double  angles  used  as  the  top  chord 
members.  Tests  were  performed  on  single  angles  and  on  double  angles 
connected  to  act  as  a  unit.  All  failures  were  associated  with  flexural 
buckling  or  local  failure,  no  torsional  failures  were  reported.  Lenzen's 
results  cannot  be  directly  compared  to  column  tests  as  they  represent 
conditions  of  loading  and  end  restraint  peculiar  to  joist  chord  segments. 
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There  is  a  lack  of  experimental  data  on  the  buckling  strength 
of  double  angle  sections.  In  addition  the  buckling  strength  of  such 
members  in  the  inelastic  range  has  not  been  determined  analytically, 
and  in  fact,  residual  stress  distributions  for  angle  sections  have  not 
been  reported. 
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Figure  2.2  Buckling  Deformations 
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Figure  2.3  Residual  Stress-Distribution-Single  Angle  (Ref.  3) 
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CHAPTER  III 

MATERIAL  PROPERTIES 

A  wide  range  of  angle  sizes  are  manufactured  with  leg  lengths 

varying  from  1"  to  9"  and  thicknesses  from  1/8"  to  1-1/8".  The  more 

common  sizes  are  generally  available  in  Canada  in  three  grades  of  steel; 

ASTM-A36 ,  CSA-G40. 1 2  ,  and  ASTM-A441 . 

Material  properties  were  obtained  from  tests  on  five  different 

angle  sections  of  G40.12  steel.  The  five  sections  listed  in  Table  3.1 

were  selected  to  provide  a  large  range  of  leg  lengths  and  to  be  large 

enough  to  provide  easy  sectioning. 

The  chemical  composition  and  mill  test  results  are  given  in 

Table  3.1.  These  results  conform  to  CSA-G40.12  specifications^.  Mill 

test  results  were  not  available  for  the  4  x  7  x  3/8  angle. 

Tables  3.2(a)  to  (c)  summarize  the  material  properties  obtained 

from  laboratory  tension  tests  performed  on  the  milled  residual  stress 

specimens.  The  location  of  the  test  specimens  in  the  section  are 

shown  in  the  inset.  In  this  table  a  denotes  the  static  yield  stress 

and  a  the  ultimate  stress;  e  represents  the  calculated  yield  strain 
u  y 

corresponding  to  a  modulus  of  elasticity  of  E  =  29,600  ksi^  and  the 

stress,  a  .  The  measured  strain  at  the  onset  of  strain  hardening  is 

denoted  by  £  .  and  E  .  is  the  strain  hardening  modulus.  The  tests  were 
J  st  st 

performed  in  a  Baldwin  hydraulic  machine  using  flat  tensile  specimens 
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similar  in  shape  to  those  required  by  ASTM  specification  A370-65  . 

The  static  yield  values  were  taken  after  a  five  minute 
period  at  zero  strain  rate.  Values  of  E  ^  were  obtained  by  graphically 
measuring  the  initial  slope  of  the  strain-hardening  portion  of  the 
load-strain  curve.  The  plastic  region  of  the  load-strain  curve  was 
pronounced  in  the  large  majority  of  the  tests  with  the  exception  of 
several  specimens  from  the  heel  of  the  angles,  which  exhibited  an 
elastic-strain  hardening  response. 

Residual  stress  distributions  were  obtained  by  the  method 

of  sectioning^.  Three  test  samples  were  cut  from  each  angle  section 

and  the  results  are  plotted  individually  in  Figures  3.1(a)  to  3.1(e). 

Plotted  stresses  were  measured  as  strains  and  converted  to  stresses 

using  E  =  29,600  ksi.  Curvature  effects  due  to  bowing  upon  sectioning 
20 

were  checked  ,  and  found  to  be  negligible.  The  stress  distributions 
for  each  angle  section  were  similar  with  excellent  agreement  for  the 
three  tests  performed  on  a  single  angle  size. 

Symmetrical  residual  stress  patterns  were  not  observed  even 
for  the  equal  legged  angles.  Maximum  compressive  stresses  occurred  at 
or  near  the  heel  of  the  angle  in  all  tests.  The  maximum  compressive 
residual  stress  measured  was  14  ksi. 

The  inset  in  Table  3.3  shows  the  assumed  residual  stress 
patterns  used  for  later  computations  of  the  buckling  strength.  The 
assumed  pattern  places  maximum  compressive  residual  stress  values  at 
the  angle  heel  and  leg  tips  and  envelopes  the  measured  values  in  most 
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cases.  The  maximum  tensile  residual  stress  is  located  at  the  mid  point 
of  each  leg. 

A  computer  program  was  developed  to  ensure  that  the  assumed 
residual  stress  distribution  did  not  produce  a  net  axial  force  or 
bending  moment  about  the  principal  axes.  The  program  accepted  the 
residual  stresses  at  the  heel  and  the  tip  of  the  short  leg  as  input 
(a  g  and  o  ,-)  and  outputed  the  remaining  peak  values. 

Table  3.3  lists  the  values  of  the  measured  (average)  residual 
stresses  at  specific  locations  and  a  reasonable  balanced  computer  solu¬ 
tion.  In  Table  3.3  compressive  stresses  are  listed  as  positive  and 
ar 2  is  taken  at  the  mid-point  of  the  shorter  leg.  In  addition,  the 
results  of  a  rather  extreme  balanced  case  (which  represents  a  maximum 
compressive  residual  stress  of  13  ksi  (as  implied  in  Reference  5)  are 
tabled.  Figure  3.2  shows  a  comparison  of  the  three  stress  distributions 
for  the  4  x  4  x  3/8"  angle. 
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Figure  3.1b  Residual  Stress  Measurements 
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Figure  3.1c  Residual  Stress  Measurements 


26 


COMP. 
(KS I ) 


i  i _ I _ 1 _ I _ l 

10  5  0  5  10  15 


TENS.  COMP. 

(KS I )  (KS I ) 


Figure  3. Id  Residual  Stress  Measurements 
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Figure  3.1e  Residual  Stress  Measurements 
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ANGLE 

RESIDUAL 

DISTRIB. 

1 

IKS!) 

<7  r2 

(KSI) 

cr  r3 

(KSI) 

4 

(KSI) 

°>5 

(KSI) 

4*  3*3/8 

MEASURED 

+  2.0 

-3.0 

+  5.0 

-5.0 

+  1.0 

BALANCED  VALUE 

+  5.7 

-5.6 

+  5.0 

-4  .2 

+  4.0 

EXTREME  BAL¬ 
ANCED  VALUE 

+13.0 

-13.0 

+13.0 

-13.0 

+13.0 

4*4*3/8 

MEASURED 

+  2.0 

-  1.0 

+  5.0 

-2.0 

+  1.0 

BALANCED  VALUE 

+  3.0 

-  3.0 

+  3.0 

-3.0 

+  3.0 

EXTREME  BAL¬ 
ANCED  VALUE 

+  13.0 

-13.0 

+  13.0 

-13.0 

+  13.0 

4*  6*  3/8 

MEASURED 

+  2.0 

-  5.0 

+  2.0 

-  1.0 

+  1.0 

BALANCED  VALUE 

+2.3 

-2.4 

+  3.0 

-3.8 

+  4.0 

EXTREME  BAL¬ 
ANCED  VALUE 

+  13.0 

-13.0 

+  13.0 

-13.0 

+13.0 

4*7*3/8 

MEASURED 

+  0.0 

-4.0 

+  10.0 

-  1.0 

+  4.0 

BALANCED  VALUE 

+  7.1 

-6.9 

+  6.0 

-4.3 

+  4.0 

EXTREME  BAL¬ 
ANCED  VALUE 

+  13.0 

-13.0 

+13.0 

-13.0 

+13.0 

4*  4*  1/4 

MEASURED 

-2.0 

-3.0 

+  6.0 

-3.0 

+3.0 

BALANCED  VALUE 

+  3.0 

-3.0 

+  3.0 

-3.0 

+  3.0 

EXTREME  BAL¬ 
ANCED  VALUE 

+  13.0 

-13.0 

+13.0 

-13.0 

+13.0 

*  COMPRESSION  POSITIVE  (+) 


TABLE  3.3  Residual  Stress  Distribution  Values 
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CHAPTER  IV 

ANALYTICAL  INVESTIGATION 


The  differential  equations  expressing  equilibrium  of  an 
axially  loaded  singly  symmetric  column  in  the  deformed  position  were 
given  as  equations  2.6  to  2.8  in  Chapter  II.  For  simply  supported 
flexural  and  torsional  boundary  conditions,  the  condition  that  the 
determinant  of  the  coefficient  matrix  vanishes  (buckling  condition) 
results  in  the  solutions  given  by  equation  2.9.  Of  the  three  possible 
solutions,  two  are  real  and  are  given  by: 


P 


x 


(4.1) 


and  P,  equal  to  the  root  of: 

< py 


-  £  (P-Py)(P-P^)  -  PyQ2  (4.2) 

The  evaluation  of  equations  4.1  and  4.2  requires  known  values 
of  the  section  properties  and  the  column  length,  to  calculate  the  two 
critical  buckling  loads,  Px  and  P^y.  In  the  elastic  range  the  critical 
loads  corresponding  to  a  given  length  can  be  determined  directly,  but 
in  the  inelastic  range  the  section  properties  depend  on  the  extent  of 
yielding  across  the  section,  which  in  turn  is  influenced  by  the  applied 
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load.  Thus  an  indirect  approach  must  be  used  to  solve  Equations  4.1  and 


4.2. 


An  elastic-plastic  stress-strain  curve  (Figure  4.1)  is  assumed 


for  the  material.  The  cross-section  is  divided  into  segments  and  a 
uniform  strain,  e  ,  is  applied  to  the  entire  cross-section.  The  total 

a 

strain  in  any  segment,  ,  is  then  equal  to  the  sum  of  the  applied 
strain,  e  ,  and  the  residual  strain,  e  . 


r 


(4.3) 


The  stress-strain  curve  can  be  entered  with  the  total  strain  and  the 
corresponding  stress  in  the  segment,  a,  is  determined  as  shown  in 
Figure  4.1.  The  total  load,  P,  on  the  cross-section  corresponding  to 
the  applied  strain  level  and  the  assumed  residual  strain  pattern  can  be 
calculated  as: 


(4.4) 


P  =  adA 
A 


The  factor  K  which  arises  from  the  differential  warping  of  two  adjacent 
cross-sections  can  be  calculated  as: 


K  =  aa2dA 
A 


(4.5) 


where  a  represents  the  distance  of  the  segment  under  consideration 
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from  the  shear  centre. 

At  a  zero  applied  strain  level  the  value  of  K  represents  the 
reduction  in  torsional  buckling  strength  due  to  the  residual  strain 
distribution.  This  reduction  is  applicable  in  both  the  elastic  and 
inelastic  ranges. 

The  flexural  and  torsional  section  properties  are  calculated 
by  assuming  that  no  resistance  is  offered  by  the  yielded  segments.  All 

o 

section  property  calculations  are  based  on  thin  wall  theory  .  The 
calculation  of  the  St.  Venant  torsional  resistance,  Ky,  is  performed 

o 

on  the  original  elastic  cross-section  .  For  a  given  level  of  applied 
strain,  e^,  the  corresponding  applied  load  and  section  properties  are 
then  known.  Equation  4.1  and  4.2  can  then  be  solved  directly  to  ob¬ 
tain  the  lengths  corresponding  to  flexural  and  lateral-torsional 
buckling,  Lx  and  ,  respectively. 

A  computer  program  listed  in  Appendix  A,  was  developed  to 
calculate  the  flexural  and  torsional  section  properties  and  solve 
the  buckling  equations  at  various  applied  strain  levels.  A  simplified 
flow  chart  for  the  computer  program  is  given  in  Figure  4.4. 

The  warping  moment  of  inertia,  I  ,  was  computed  by  assuming 
that  the  two  angles  were  joined  on  the  gauge  line  by  a  thin  strip,  as 

9 

shown  in  Figure  4.3.  The  value  of  I  for  a  single  angle  is  zero  thus 
the  angle  sections  must  act  as  a  unit  for  the  warping  moment  of  inertia 
to  have  a  finite  value. 
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UNIFORM  APPLIED  STRAIN 
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Figure  4.2 


Strain  Patterns 
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Figure  4.3  Double  Angle  Connection  at  Gauge  Location 
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Flow  Chart  for  Computer  Program 
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CHAPTER  V 

TEST  PROGRAM  AND  RESULTS 

The  analysis  described  in  the  previous  chapters  proceeds  by 
first  selecting  the  level  of  applied  strain  to  be  considered.  The 
applied  strain  is  superimposed  on  the  residual  strain  in  each  segment 
of  the  cross-section  and  the  corresponding  stress  is  determined.  The 
section  properties  are  then  computed  using  a  numerical  integration 
procedure;  the  yielded  segments  are  omitted  from  this  calculation. 
Approximately  eighty  segments  per  leg  were  used  to  obtain  an  accurate 
solution.  The  stresses  on  each  segment  are  then  summed  to  determine 
the  applied  load  and  Equations  4.1  and  4.2  are  solved  for  the  corres¬ 
ponding  equilibrium  lengths.  In  addition,  Equation  2.10  is  solved  for 
the  length  corresponding  to  (fictitious)  y-y  axis  buckling. 

This  procedure  is  repeated  for  increasing  values  of  the  ap¬ 
plied  strain  level,  corresponding  to  higher  buckling  loads  and  pro¬ 
gressive  yielding  of  the  member.  Finally,  the  complete  cross-section 
has  yielded  and  the  corresponding  buckled  lengths  are  equal  to  zero. 

The  results  of  the  analyses  are  plotted  in  the  form  of 
column  curves.  In  these  curves  the  stress  corresponding  to  buckling, 
a,  is  non-dimensionalized  as  a /a  ,  where  represents  the  yield 
stress,  and  is  plotted  against  the  slenderness  ratio,  L/r  .  The  value 
of  a  is  given,  at  each  strain  level,  by: 
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a  =  Z 


aAA 

A 


(5.1) 


where  a  represents  the  stress  on  a  particular  segment  of  the  cross- 
section  and  AA  represents  the  area  of  the  segment. 

The  results  are  presented  in  Figures  5.1  to  5.5.  In  each 
case  two  curves  are  given;  that  corresponding  to  x-x  axis  flexural 
buckling  and  that  corresponding  to  lateral-torsional  buckling. 

Each  curve  has  two  portions;  the  elastic  (concave)  portion  pertaining 
to  large  values  of  L/r  and  the  inelastic  (convex)  portion  for  smaller 

A 

slenderness  ratios.  In  general,  the  curves  corresponding  to  flexural 
buckling  are  shown  as  dashed  and  those  corresponding  to  lateral-torsional 
buckling  as  solid  curves. 

The  curves  corresponding  to  (fictitious)  y-y  axis  buckling 

are  not  shown.  In  all  cases,  however,  they  lie  just  above  the  curves 

corresponding  to  1 ateral -torsional  buckling.  For  a  given  slenderness 

ratio,  the  critical  stresses  are  within  2%  of  one  another.  Those 

shown  in  Table  5.1  coincide  in  the  second  decimal  place. 

Influence  of  Angle  Symmetry  -  Figures  5.1(a)  through  5.1(e)  show  column 

curves  for  five  double  angle  sections.  The  sections  are  composed  of 

115  3 

angles  ranging  from  the  3^  x  2^-  x  to  a  7  x  4  x  g".  In  addition  to 

variations  in  size  the  angles  have  differing  ratios  of  r  to  r  ;  where 

x  y 


r  denotes  the  radius  of  gyration  of  the  section  about  the  x  axis  and 

X 


3  ii 


r  the  radius  of  gyration  about  the  y  axis.  For  example,  the  7  x  4  x  g 1 

3 

angle  has  a  ratio  of  1.43  while  the  4  x  4  x  g "  has  a  ratio  of  0.66.  Each 
of  the  curves  shown  in  Figure  5.1  was  obtained  for  a  member  having  a 
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yield  stress  level  of  44  ksi,  a  maximum  compressive  residual  stress  of 
13  ksi  and  the  section  is  formed  by  placing  the  long  legs  of  the  angles 
back  to  back.  The  separati on, assumed  to  be  caused  by  a  gusset  plate, 
was  related  to  the  leg  thickness.  In  the  analysis,  the  angles  were 
assumed  to  act  as  a  unit,  connected  at  the  gauge  line  (Table  1.1). 

The  members  having  large  values  of  r  (r  /r  ratios  £  1.0) 

y  ^  y 

have  lateral-torsional  buckling  strengths  greater  than  flexural  buck¬ 
ling  strengths.  As  r  /r  increases  above  unity  the  trend  reverses. 

x  y 

This  is  shown  for  example,  by  a  comparison  of  the  curves  for  the 
4  x  4  x  q- "  (r  /r  =  0.66)  angle  (Figure  5.1(c))  with  curves  for  the 
7  x  4  x  I""  (r  /r  =  1.43)  angle  (Figure  5.1(e)).  This  is  as  expected 

sj 

since  the  lateral-torsional  buckling  loads  and  the  y-y  flexural  buckling 
loads  are  virtually  coincident.  Note  that  in  some  cases  cross-overs 
occur  in  the  inelastic  range  (for  example  see  Figure  5.1(a)). 

Influence  of  Leg  Thickness  -  Figures  5.2(a)  through  5.2(e)  show  column 
curves  for  double  angle  sections.  The  sections  are  composed  of  6"  x  4" 
angles  ranging  in  thickness  from  5/16"  to  3/4".  The  radius  of  gyration 
ratios  fall  within  a  relatively  narrow  range,  from  1.20  for  the  5/16" 
leg  thickness  to  1.05  for  the  3/4"  thick  angle  section.  The  angles 
are  placed  with  the  long  legs  back  to  back  and  the  assumed  gusset 
plate  thickness  is  compatible  with  the  leg  thickness  of  the  angles. 

The  members  analyzed  had  a  yield  stress  of  44  ksi  and  a  maximum  com¬ 
pressive  residual  stress  of  13  ksi. 

The  members  all  have  small  values  of  ry(rx/ry  ratios  >1.0) 
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which  produces  flexural  buckling  strengths  greater  than  lateral- 
torsional  buckling  strengths.  Note  that  in  some  cases  cross-overs 
do  occur  in  the  inelastic  range  (Figure  5.2(d)). 

Influence  of  Angle  Separation  -  Figure  5.3  shows  the  column  curves  for 
double  angle  sections  composed  of  6  x  4  x  g"  angles  having  the  long 
legs  back  to  back.  The  gusset  plate  thickness  (separation  of  the  long 
legs)  ranges  from  3/8"  to  3/4".  The  radius  of  gyration  ratios  range 
from  1.19  for  the  3/8"  separation  to  1.10  for  the  3/4"  separation. 

The  yield  stress  level  and  residual  stress  distribution  were  the 
same  as  the  previous  analyses.  The  flexural  buckling  curves  for 
the  four  members  coincide.  The  lateral-torsional  strengths  vary, 
however,  as  shown  by  the  solid  curves. 

The  members  again  illustrate  dominant  flexural  buckling 
strengths  with  lateral-torsional  buckling  strengths  approaching  the 
flexural  strengths  as  the  ratio  rx/r  approaches  unity. 

Influence  of  Yield  Stress  Level  -  Figures  5.4(a)  through  5.4(d)  show 

3 

column  curves  for  double  angle  sections  composed  of  6  x  4  x  g"  angles 
having  the  long  legs  back  to  back,  and  a  constant  gusset  plate  thick¬ 
ness  of  4".  The  yield  stress  levels  were  varied  to  represent  four 
common  materials;  36  ksi  yield-ASTM-A36 ,  44  ksi  yield-CSA-G40.12,  50 
ksi  yiel d-ASTM-A441  and  100  ksi  yield-ASTM-A514.  The  maximum  com¬ 
pressive  residual  stress  was  assumed  to  be  13  ksi  for  all  members. 

The  flexural  buckling  strengths  were  greater  than  the 
lateral-torsional  strengths  in  all  cases  as  the  ratio  rx/ry  was 


■ 


44 


greater  than  unity.  An  increase  in  the  yield  stress  level  effectively 
increases  the  elastic  limit,  since  the  ratio  of  maximum  compressive 
residual  stress  to  yield  stress  is  reduced.  For  example  ASTM-A36  steel 
has  an  elastic  limit  of  0.64  a  corresponding  to  a  slenderness  ratio 

V 

of  approximately  100  for  lateral-torsional  buckling,  whereas  ASTM-A514 
has  an  elastic  limit  of  0.87  a  corresponding  to  a  slenderness  ratio 
of  approximately  50. 

Influence  of  Residual  Stress  Distribution  -  Figure  5.5  shows  the  column 

3 

curves  for  double  angle  sections  composed  of  7  x  4  x  g"  angles  having 

5 

the  long  legs  back  to  back  and  a  constant  gusset  plate  thickness  of  g-'.' 

The  yield  stress  level  was  held  constant  at  44  ksi ,  with  three  different 

3 

residual  stress  distributions.  The  7  x  4  x  g  angle  was  used  as  this 
represented  the  section  with  the  maximum  measured  residual  stress  (14 
ksi  compressive).  The  first  distribution  corresponded  to  the  unbalanced 
measured  values;  these  were  then  balanced  by  a  computer  solution 
(Chapter  IV)  for  the  second  distribution.  The  third  balanced  distri¬ 
bution  was  based  on  the  recommendations  of  CSA-S16  and  has  a  maximum 
compressive  residual  stress  of  13  ksi. 

The  flexural  buckling  strengths  were  greater  than  lateral- 
torsional  buckling  strengths  for  these  members,  as  the  ratio  of  rx/ry 
exceeded  unity.  The  unbalanced  measured  residual  stress  distribution 
produced  the  maximum  buckling  strengths  for  both  buckling  modes,  although 
the  curves  showed  a  sudden  break  when  a  substantial  portion  of  the  cross- 
section  had  yielded.  The  balanced  measured  distribution  smoothed  out 
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the  column  curves  and  the  ultimate  strengths  were  below  those  ob¬ 
tained  with  the  unbalanced  distribution.  The  13  ksi  maximum  com¬ 
pressive  residual  stress  distribution  was  below  the  other  corresponding 
curves  and  produced  a  conservative  estimate  of  the  buckling  strengths. 
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Figure  5.1a  Influence  of  Member  Asymmetry 


4  x  3  x  3/8"  ANGLES 
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Figure  5.1b  Influence  of  Member  Asymmetry 
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Figure  5.1c  Influence  of  Member  Asymmetry 


6x4  x  3/8"  ANGLES 
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Figure  5. Id  Influence  of  Member  Asymmetry 
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Figure  5.1e  Influence  of  Member  Asymmetry 
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Figure  5.2a  Influence  of  Leg  Thickness 
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Figure  5.2b  Influence  of  Leg  Thickness 


6  x  4  x  i/2"  ANGLES 
crv  =  44  ksi  “1  r 
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Figure  5.2c  Influence  of  Leg  Thickness 
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Figure  5. 2d  Influence  of  Leg  Thickness 
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Figure  5.2e  Influence  of  Leg  Thickness 
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Figure  5.3  Influence  of  Leg  Separation 
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Figure  5.4a  Influence  of  Yield  Stress 


6x  4x3/8"  ANGLES 
try  =  44 ksi  CSA -  G40.12 
<Xr  max  13  ksi 
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Figure  5.4b  Influence  of  Yield  Stress 


6  x  4  «  3/8"  ANGLES 
crv  =  50  ksi  ASTM  -  A  441 
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Figure  5.4c  Influence  of  Yield  Stress 


6  x  4  x  3/8"  ANGLES 
crv  =  100  k si  ASTM  -A514 
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Figure  5.4d  Influence  of  Yield  Stress 
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Figure  5.5  Influence  of  Residual  Stress 
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6*4* 3/8 "  ANGLE 

5/8"  GUSSET  PLATE 

LONG  LEGS  BACK  TO  BACK 

L/rx 

cr/Oy 

(x-x) 

cr/ary 

( <£y ) 

CJ/Oy 

(y-y) 

Fictitious 

25 

0.99 

0.99 

0.99 

50 

0.92 

0.92 

0.92 

75 

0.80 

0.78 

0.78 

100 

0.67 

0.50 

0.50 

1  25 

0.44 

0.34 

0.34 

1  50 

0.30 

0.23 

0.23 

1  75 

0.21 

0.17 

0.17 

200 

0.17 

0.13 

0.13 

Table  5.1  Critical  Stresses 
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CHAPTER  VI 

APPLICATION  TO  DESIGN 

Procedures  for  the  design  of  double  angle  struts  are  based 
on  the  assumption  that  flexural  buckling  will  mark  the  ultimate  strength 
of  the  member.  The  possibility  of  coupled  flexural -torsional  buckling 
is  ignored.  *  ' 

The  results  of  this  investigation  show  that  the  lateral- 
torsional  buckling  load  for  a  double  angle  strut,  P^  ,  is  less  than  the 
y-y  axis  flexural  buckling  load,  P  ,  but  the  difference  is  within  2% 
for  all  cases  investigated. 

The  virtual  coincidence  of  P .  ,  and  f>  can  be  explained  using 

4>y  y 

Equation  4.2  as  the  basis;' 

-  jr  (P-PyHP-P^,)  -  (Py0)2  -  o  (6.1) 


By  neglecting  the  differential  warping  caused  by  the  residual  stresses 

and  assuming  that  the  stress  on  a  particular  fibre  o  =  P/A,  a  quantity 
_  9 

rQ  can  be  degined  as: 


K 

P 


*o2  + 


(6.2) 


Using  this  definition,  Equation  6.1  can  be  rearranged  as: 


* 
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wi  th ; 


1 


and: 


1  *2EI 

<P  -T  [GKT  +  ~2~] 
0  L 


p,  = 


a  shape  factor,  J,  can  be  defined  as: 


0  =  [1  - 


(6.3) 


(6.4) 


(6.5) 


(6.6) 


where  J  is  positive  and  always  less  than  or  equal  to  unity. 

In  the  cases  considered  the  warping  resistance  (~~yW)  was 
significantly  less  than  the  St.  Venant  torsional  resistance  (GK^). 
Neglecting  the  warping  contribution  and  substituting  the  results  of 
Equations  6.4  to  6.6  into  Equation  6.3  the  lateral  torsional  buckling 
load  can  be  expressed  as: 


1  -  (pR)2/GKT] 

1  -  J(p?02/gkt 


(6.7) 


In  Equation  (6.7)  the  factor  J  reduces  the  value  of  the  lateral-torsional 
buckling  load  below  the  (fictitious)  y-y  axis  flexural  buckling  load. 
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For  example,  consider  a  6  x  4  x  |"  double  angle  strut  with 
1/2  gusset  plate  and  the  long  legs  placed  back  to  back.  For  a  slender¬ 
ness  ratio,  L/rx,  of  100  the  warping  torsional  resistance  is  approxi¬ 
mately  4%  of  the  St.  Venant  resistance  and  can  be  neglected.  The 

corresponding  value  of  J  is  0.97;  this  implies  that  P,  will  be 

<\>y 

approximately  2%  less  than  P  .  In  the  inelastic  range  at  a  slender¬ 
ness  ratio,  L/rx,  of  50  the  difference  between  the  two  loads  is  less 
than  1%.  These  results  are  typical  of  the  cases  considered. 

In  the  design  of  a  double  angle  strut,  the  column  strength 
will  be  based  on  the  lesser  of  the  y-y  and  x-x  axes  buckling 
strengths.  Since  the  y-y  axis  buckling  strength  virtually  coincides 
with  the  lateral-torsional  buckling  strength,  this  procedure  results  in 
a  correct  assessment  of  the  member  strength. 

The  allowable  stress  procedure  for  column  design  restricts 
the  computed  stress  on  the  member  (under  design  loads)  to  a  value  ob¬ 
tained  by  dividing  the  critical  stress  (member  strength)  by  a  factor 
of  safety.  The  factors  of  safety  used  in  North  American  building 
codes  are  based  on  the  buckling  strength  of  rolled  wide-flange  members 
and  have  been  checked  against  large-scale  tests  on  such  members. 

The  residual  stress  distribution  for  a  double  angle  strut 
and  its  influence  on  the  gradual  plastification  of  the  section  may  not 
be  the  same  as  that  for  a  rolled  wide-flange  section.  In  order  to 
assess  the  suitability  of  the  above  design  provisions  for  double  angle 
struts,  the  buckling  strengths  obtained  from  the  computer  analysis  are 
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compared  with  the  provisions  of  the  CSA^  and  AISC^  standards  in 
Figures  6.1  and  6.2. 

Figure  6.1  and  6.2  plot  the  flexural  and  lateral-torsional 
buckling  strengths  respectively,  for  a  4  x  3  x  double  angle  strut, 
placed  with  long  legs  back  to  back  and  separated  by  a  1/2"  gusset 
plate.  The  material  is  CSA-G40.12  steel  with  a  specified  minimum 
13  ksi  (using  the  assumed  distribution  of  Chapter  III). 

Figure  6.1  plots  the  critical  stress  for  x-x  axis  flexural 

buckling,  a/a  ,  against  the  slenderness  ratio,  L/r  ;  whereas  the 
y  x 

curves  representing  1 ateral -torsional  buckling.  Figure  6.2,  plot 
the  critical  stress  against  the  y-y  axis  slenderness  ratio,  L/r  ,  to 
simplify  the  comparisons.  The  computer  solutions  are  shown  as  heavy 
solid  curves  with  the  CSA  allowable  stresses  as  lighter  solid  curves 
and  the  AISC  allowable  stresses  as  light  broken  lines.  The  factors  of 
safety,  based  on  comparisons  of  the  flexural  and  lateral  buckling 
strengths  with  the  CSA  and  AISC  design  provisions  are  shown  in  Table 
6.1 . 

In  the  elastic  range  these  factors  of  safety  are  constant  at 
1.96.  This  value  is  normally  taken  as  1.92;  the  difference  is  the 
result  of  using  E  =  29,600  ksi  in  the  computer  program  instead  of  the 
more  conservative  value  of  E  =  29,000  ksi  *  . 

In  the  inelastic  range  the  factors  of  safety  provided  by  the 
CSA  provisions  decrease  gradually  from  1.96  at  the  elastic  limit  to 
1.67  at  a  slenderness  ratio  of  zero.  The  factors  of  safety  provided 
by  the  AISC  provisions  increase  from  1.96  at  the  elastic  limit  to 
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2.06  at  a  slenderness  ratio  of  100,  then  decrease  gradually  to  1.67 
at  a  slenderness  ratio  of  zero. 

In  Table  6.1,  the  critical  stresses,  a  ,  as  given  by  the 

C  i 

20 

Column  Research  Council's  ,  basic  column  formula,  are  also  listed; 
where 

acr  *  V1  -  ^  (1  ‘  ^)(^)(T>2]  (6-8) 

Using  a  peak  compressive  residual  stress,  c  of  13  ksi ,  the  predictions 
of  the  CRC  equation  matched  the  computed  flexural  and  lateral-torsional 
buckling  stresses  to  within  1%. 

Local  buckling  of  the  angle  was  not  considered  in  this  in¬ 
vestigation.  Column  curves  have  been  presented  for  angle  sections 
having  leg  width-to-thickness  ratios,  b/t,  exceeding  the  CSA  require¬ 
ment  of  75//F  and  the  AISC  requirement  of  76//F  ,  in  order  to  investi- 
gate  the  significance  of  the  various  structural  parameters.  Where  the 
leg  width-to-thickness  ratio  exceeds  the  code  limitation  a  reduction 
factor  would  be  applied  to  the  allowable  stress  to  account  for  the 
possibility  that  local  buckling  of  the  plates  which  make  up  the  cross- 
section  may  occur  before  overall  buckling  of  the  member. 
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Figure  6.1  Ultimate  and  Design  Curves  Flexural  Buckling 


4  x  3  x  3/8"  ANGLES 
cty  =44  ksi 

o*r  max  =  1 3  ksi 
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Figure  6.2  Ultimate  and  Design  Curves  Lateral-Torsional  Buckling 
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L/rx 

COMPUTED  1 
ULTIMATE 
STRENGTH 
(KSI) 

COLUMN 

RESEARCH 

COUNCIL 

(KSI) 

CSA-S16 

(DESIGN) 

(KSI) 

AISC 

(DESIGN) 

(KSI) 

COMPUTED  FACTOR 
OF  SAFETY 

ULT,  STR. 

ULT.  ST.R. 

CSA-S16 

AISC 

0 

44.0 

44.0 

26.4 

26.4 

1.67 

1.67 

25 

43.5 

43.2 

25.7 

24.5 

1.69 

1.78 

50 

40.0 

40.0 

22.3 

21.8 

1.80 

1.83 

75 

35.6 

35.6 

19.0 

18.4 

1.88 

1.94 

100 

29.2 

29.2 

14.9 

14.2 

1.96 

2.06 

125 

18.8 

18.8 

9.6 

9.6 

1.96 

1.96 

150 

12.9 

12.9 

6.6 

6.6 

1.96 

1.96 

175 

9.6 

9.6 

4.9 

4.9 

1.96 

1.96 

L/ry 

ii 

ii 

ii 

n 

it 

ii 

0 

44.0 

44.0 

26.4 

26.4 

1.67 

1.67 

25 

43.7 

43.2 

25.7 

24.5 

1.70 

1.78 

50 

41.8 

40.5 

22.3 

21.8 

1.87 

1.92 

75 

36.3 

36-4 

19.0 

18.4 

1.91 

1.97 

100 

29.2 

29.2 

14.9 

14.2 

1.96 

2.06 

125 

18.8 

18.8 

9.6 

9.6 

1.96 

1.96 

150 

12.9 

12.9 

6.6 

6.6 

1.96 

1.96 

175 

9.6 

9.6 

4.9 

4.9 

1.96 

1.96 

Table  6.1  Factors  of  Safety 
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CHAPTER  VII 

SUMMARY  AND  CONCLUSIONS 

This  investigation  attempts  to  define  the  flexural  and 
1  ateral -torsional  buckling  strengths  of  double  angle  struts. 

Residual  stress  distributions  were  measured  on  five  angle 
sections  of  CSA-G40.12  steel  and  are  presented  in  Figures  3.1(a)  to 
3.1(3);  along  with  the  corresponding  material  properties. 

The  residual  stress  distribution  was  idealized  and  used 
in  a  computer  program  to  predict  the  flexural  and  lateral-torsional 
buckling  strengths.  The  cross-section  properties  calculated  were 
based  on  the  elastic  core  of  the  cross-section  remaining  at  each  stage 
of  loading.  In  the  computations,  the  two  angles  composing  the  member 
were  assumed  to  be  joined  by  a  thin  strip  of  material  located  at  the 
gauge  line.  This  implies  that  the  two  angles  act  as  a  unit  to  resist 
the  buckling  motions. 

Column  curves  are  plotted  as  Figures  5.1  to  5.5  and  illu¬ 
strate  the  influence  of;  member  asymmetry,  leg  thickness,  leg  separation, 
yield  stress  and  residual  stresses  on  the  flexural  and  lateral-torsional 
buckling  strengths  of  the  member. 

In  all  cases  investigated  the  lateral-torsional  buckling  load, 
?±  ,  was  within  2%  of  the  y-y  axis  flexural  buckling  load.  This  is 

<py 

because  the  shape  factor,  J,  for  double  angle  sections  is  approximately 


-erpnfc  avii  9rtt  ,®no* 
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equal  to  one,  throughout  the  complete  loading  range. 

The  flexural  and  lateral-torsional  buckling  strengths  compare 
favorably  with  predictions  given  by  the  CRC  basic  column  strength 
equation.  In  addition  the  same  factors  of  safety  required  by  the  CSA 
and  AISC  building  codes  for  rolled  wide  flange  sections  are  achieved 
by  double  angle  sections  designed  only  on  the  basis  of  flexural 
buckl i ng. 

The  design  of  double  angle  struts  considering  only  the  strong 
and  weak  axis  buckling  strengths  is  satisfactory,  since  the  lateral- 
torsional  buckling  strength  is  within  2%  of  the  y-y  axis  buckling 
strength  normally  computed  by  the  designer. 


' 


NOMENCLATURE 


Cross  sectional  area 
Modulus  of  elasticity 
Strain-hardening  modulus 
Factor  of  safety 
Gauge  location 
Shear  modulus 
Moment  of  inertia 
Warping  moment  of  inertia 
Shape  factor 

Differential  warping  reduction 
St.  Venant  torsional  constant 
Length 

Applied  load 
Critical  buckling  load 
Flexural  buckling  load  X-X  axis 
Flexural  buckling  load  Y-Y  axis 
Lateral-torsional  buckling  load 
Torsional  buckling  load 
Radius  of  gyration 
Polar  radius  of  gyration 
Displacement  in  x  direction 


if  I  '■ 


>ol  H  ~ik  sn  s 
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v  Displacement  in  y  direction 

Xq  X  coordinate  of  shear  centre 

Yq  Y  coordinate  of  shear  centre 

e  Strain 

e  Applied  strain 

3 

e  Residual  strain 

r 

e  ^  Strain  at  onset  strain-hardening 

Total  strain 
e  Yield  strain 

y 

<j)  Rotati  on 

a  Stress 

a  Applied  stress 

3 

a  Critical  buckling  stress 

c 

a  Residual  stress 

r 

ct  Total  stress 

a  Yield  stress 

y 


from  centroid 
from  centroid 


, 

-» 
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A  -  1 


LEVEL  17  (  1  NOV  68)  OS/360  FORTRAN  H 

COMPILER  OPTIONS  -  NAME=  MA I N , 0PT=02 , L I NECNT  =  5 9 , SOURC E , EBC 0 1 C , NQL I S T , NOOECK , L OAD^l 
INTEGER  OUTPUT  ,  D1  ,  02,  03,  M  ,  N,  H  NO  MAP , NOED l T ,  ID , NOXRE F| 

C  DEFINE  CARD  READER  NUMBER 

INPUT  =  5 

C  DEFINE  PRINTER  UNIT  NUMBER 

OUTPUT  =  6 

PEAL  L 1 , L2 , L3 , LOAD , I  XX , I  YY  ,  LRX ,  L  RY  ,  I  WY  ,  I WX  ,  I W  ,  L LI ,  LL2 

DIMENSION  SR1I90)*SR2(90)»SR3I5),WO(500)*WN(500),WNN(500)*AIW(500) 

LOAD  =  0.0 
IXX  =  0.0 
I  YY  =  0.0 
A  =  0.0 
AD  =  0.0 

C  READ  MATERIAL  PROPERTIES  IN  KSI. 

READ! INPUT, 1 )E,STPESY,G 

1  FORMAT ( 3F 10 . 2 ) 

SY  =  STRESY  /  E 

C  READ  MEMBER  SIZES  IN  IN.  AND  DIVISIONS  REQUIRED  (AS  PER  CSA-S16) 

READ ( INPUT ,2 )  L1,T1,D1,L2,T2,D2,L3,T3,D3,K 

2  FORMAT  (  3  I2F5.2.I5)  ,12) 

C  REDEFINE  NEW  LENGTHS  ADJUST  DIMENSIONS  TO  CENTRE  LINE 
Ll=Ll-T2/2. 

L2=L2-Tl/2. 

L  3=L  3+T  2 
AL 1=L 1/DI 
AL  2  =  L2/02 
AL  3  =  L3/D3 
GAUGE=K* AL2 

C  READ  ALL  RESIDUAL  STRAINS  INDIVIDUALLY 

READ( INPUT, 3  )  ( SR  1 ( I ) , I  =  1 , D 1 ) 

3  FORMAT  (  8F10.8  ) 

READ( INPUT ,3 ) ( SR2( I) ,I=1,D2) 

RE AD ( INPUT, 3  )  ( SR 3 (  I ) , I = 1 , D3 ) 

C  PRINT  HEADINGS  AND  DATA  INPUT 

C  READ  ONE  VALUE  OF  APPLIED  STRAIN 

10  READ! INPUT, 4, END=999)SA 

4  FORMAT  (  F10.8) 

WRITEIOUTPUT ,20) 

20  FORM  AT ( •  1*  , 51X ,  'CIVIL  ENGINEERING  DEPARTMENT'  /  54X  , 

l'UNIVERSITY  OF  ALBERTA'  //  32X  , 'ANALYSIS  OF  RESIDUAL  STRAINS  IN  D 
20UBLE  ANGLE  COMPRESSION  MEMBERS'  /) 

WRITEIOUTPUT, 21  I  E  ,  STRESY  ,  SY 

21  FORMAT!  9X,  'MATERIAL  PROPERTIES'  //  14X  ,  'MODULUS  OF  ELASTICITY 

1  =  ',  F10.2  ,  '  KSI.'  //  14X  ,  'YIELD  STRESS',  11X,'=  '  ,F10.2  , 

2  'KSI.'  //  1 4X  , 'YIELD  STRAIN',  liX,  •=  ',  F10.8  ,'  IN/IN'  /) 
WRITEIOUTPUT  ,2  2)  0 1 , T  1  , L 1 , D2 , T2 , L 2 , D3 , T 3 , L 3 , GAUGE 

22  FORMAT!  9X,  'STRUCTURF  TOPOLOGY*  //  20X  , 'MEMBER',  18X,  'NUMBER  OF 
1'  ,  12X  , 'THICKNESS'  ,12X  ,  'LENGTH'  /  16X  ,  'NAME'  ,6X  , 'NUMBER', 


2 

12X,  'DIVISIONS 

'  ,  14X  , 

•INCHES' 

• 

1 3  X  , 

•  INCHES 

•  // 

1  4  X 

, 

3 

•FREE  LEGS 

1'  , 1 8X  , 

I  3 

,  17X 

t 

F  5 . 3 

,  1  4X  , 

F  5 . 3 

// 

14X  , 

4 

•JOINED  LEGS 

2'  ,  1 8  X  , 

13 

,  17X 

t 

F  5 . 3 

,  1 4X  , 

F  5 . 3 

// 

1 4  X , 

5 

•  JOINER 

3*  , 18X  , 

13 

,  17X 

f 

F  5 . 3 

,  14X  , 

F  5.  3 

// 

14X, 

6 

•GAUGE  LENGTH  = 

',  F5.3  , 

i 

IN' 

// 

) 

C  GAUGF  LENGTH  MEASURED  FROM  CENTRE  LINE  SECTION  ONE 
WRITEIOUTPUT, 23) 

23  FORMAT  (  9X,  'RESIDUAL  STRAINS  (COMPRESSION  POSITIVE)'//  14X  , 
1  51 'MEM  DIV  STRAIN  •  )  ) 

WRITEIOUTPUT  ,25)1  I  I,  SR  1  I  I  )  )  ,  I  =  1 ,  D 1 ) 

25  FORMAT  (15X  ,  5('l',  2X,I3,2X,  F10.8,  4X  )/ 

1(15X,5( •  1'  ,2X, I3,2X,F10.8,4X) ) ) 

WRITE (OUTPUT,2  7) I  I  I  , SR2 II) ) , I  =  1 , D2 ) 

27  FORMAT  ( 15X  ,  5('2‘  ,2X,I3  ,  2X,  F1C.8,  4X  )/ 

1 (15X,5( '2' ,2X, I3,2X,F10.8,4X) ) ) 

WRITEIOUTPUT, 29) I  I  I,  SR  3  I  I  )  )  ,  1  =  1, D3) 

29  FORMAT  (15X  ,  5('3',  2X  ,13  ,2X,  F10.8.4X  )/ 

II 15X,5( '2' ,2X, I3,2X,F10.8,4X ) ) ) 

WRI TEIOUTPUT  ,30)  SA 
WRITEIOUTPUT, 20) 

30  FORMAT!  9X  , 'APPLIED  STRAIN  =•  ,  F10.8  ,'IN/IN'  /) 

WRITEIOUTPUT, 31  ) 

31  FORMAT!  9X  , 'PORTIONS  OF  SECTION  YIELDED  DUE  TO  RESIDUAL  STRAIN' 
1  //  2 OX  , 'MEMBER'  ,  10X  , 'DIVISION*  ,  10X,  'TOTAL'  /  20X, 

2 'NUMBER '  ,  1 1 X  ,  'NUMBER'  ,  10X  ,  'STRAIN'  /  ) 

C  TEST  FOR  YIELDED  SECTIONS  AND  COMPUTE  LOCATION  OF  X-X  AXIS 

C  CALCULATE  I YY 
C  SECTION  1 

YB=L  3/2 . 

DO  50  I =1,01 
ST  =  SA  +  SR  1 1  I  ) 

IF(ST.GE.SY)ST=SY 
L0AD=L0A0*2.*ST*E*T1*AL1 
I F I  ST . GE .SY ) GO  TO  50 
XRAR=YB*L1-I I  *  AL  l l+ALl/2. 

IYY=  IYY+2 . *1  I  T  1*  A  L 1  **3/12.  )*(T1*AL1*XBAR**2)  ) 

A=A+2.*T1*AL1 

50  CONTINUE 


-  ;•  ,i' 


A  -  2 


C  SECTION  2 

DO  150  1=1,02 
ST=SA+SR2( I ) 

IF{ST.GE.SY)ST=SY 
LOAD=LOAD+2.*ST*E*T2*AL2 
I F ( ST.GE.SYIGO  TO  150 
XBAR=YB 

IYY=IYY+2.*I <AL2*T2**3  / 1 2  .  )  ♦  ( AL2 *T2*X B AR ** 2 ) ) 
A0=AD*2.*( AL2*T2*IL2-I*AL2+AL2/2. ) ) 

A=A+2.*AL2*T2 
150  CONTINUE 
C  SECTION  3 

D0250  1=1, D3 
ST=SA*SR3< I » 

IF(ST.GE.SY)ST=SY 
LOAO=LOAD4-ST*E*T3*AL3 
IF( ST.GE.SYIGO  TO  250 
XBAR=ABS( YB-I*AL3+AL3/2. ) 

I YY= I YY+ ( ( T3*AL  3**3/ 12 . )  +  ( AL 3*T3*X BAR ** 2  I  ) 
AD=AD+AL3*T3*GAUGE 
A=A+ AL3*T3 
250  CONTINUE 
YBAR= AD/A 
YAR  E A=A 

C  CALCULATE  IXX 
C  SECTION  1 

DO  300  1=1, D1 
T=T  1 

ST=SA*SR1! I ) 

I  F ( ST.GE.SY)T=0. 

D=  Y  B  AR 

IXX=IXX+2.*I  (AL1*T**3/12.  )+(ALl*T*D**2)  ) 

249  F0RMAT(5X,*T=',F12.4,3X,'D=',F12.4,3X,'IXX=',F12.4/I 
300  CONTINUE 
T  =  T2 

C  SECTION  2 

DO 3 10  1=1,02 
T=T2 

ST=S A+SR2 ( I ) 

I F ( ST.GE.SY)T=0. 

n=ABS( L2-YBAR-I *AL2+AL2/2. ) 

I  XX  =  I X  X  +  2 • *! IT*AL2**3/12. )  + ( AL 2*T*D* *2 ) ) 

310  CONTINUE 
T=T3 

C  SECTION  3 

D0320  1=1, D3 
T  =  T3 

ST=  S  A*SR  3 (  I  ) 

I E ( ST.GE.SY) T=0.0 
D=ABS(GAUGE-YBAR) 

IXX=IXX+( ( AL3*T**3/12. )+( AL3*T*D**2 ) I 
320  CONTINUE 

C  CALCULATE  FLEXURAL  BUCKLING  X  AND  Y  AXES 
RXX  =!  IXX/A)**0.5 
RYY  =< IYY/A)**0.5 

XLCR  =  ( 3. 141 59*3. 141 5 9*E* I  XX /LOAD) **0.5 
YLCR=( 3.  14159*3. 1415  9*E* I YY /LO AD ) **0 . 5 
LRX  =  XLCR  /  RXX 
LRY  =  YLCR  /  RYY 
C  OUTPUT  SECTION  PROPERTIES 

WRITEIOUTPUT ,400) 

4C0  FORMAT ( 9X  , 'PROPERTIES  OF  SECTION  NOT  YIELDED'  /) 
WRITEIOUTPUT  ,401 )  IXX  ,  I YY  ,  LOAD  ,  YBAR, 

1  RXX  ,  RYY  ,  XLCR  ,  YLCR  ,  LRX  ,  LRY 
401  FORMAT!  14X  ,'IXX  =  •  ,F7.2  ,'  IN-4'  //  14X,'IYY 

1'  IN-4*  //  14X  , 'LOAD  =  •  ,F10.2  ,'  KIPS'  //  14X 

2,  F5.2  ,  '  IN.'  // 

3 1 4X , ' R  X  =  • , F7.2  , •  IN. •  // 

414X, ' RY  =  ' , F7 . 2  , '  IN. •  // 

514X, 'CR  IT  ICAL  LENGTH  XX  AXIS  =  '  ,F7.2,'  IN'  // 

614X, 'CR I TICAL  LENGTH  YY  AXIS  =  •  ,F7.2,'  IN'  // 

7 14X , ' L/R  X  AXIS  =  '  , F  7 . 2  // 

814X , • L/R  Y  AXIS  =  ',F7.2  //  ) 

C  CALCULATE  KT  ON  ELASTIC  SECTION  ONLY 
AKT  =  0 .0 

AKT=AKTM2./3.*L1  »*(T1**3) 

AKT=AKT  +  ( 2. /3.*L2)*<  T2**3) 

AKT=AKT+( l./3.*L3)*( T3**3> 

WRITEIOUTPUT, 603IAKT 
603  FORMAT! 10X, • KT  =  »,1F10.6//) 

GUAGE  =  GAUGE 
XAXIS  =  YBAR 
OMEGA  =  0. 

IWX  =  0. 

IWY  =  0. 

B ARK  =  0 . 0 


=  '  ,  F  7. 2 , 
,'YBAR  = 


' 


o  o 


A  -  3 


C  CALCULATE  OMEGA  SECTION  1 
DO  533  I  =1,01 
OM 1  =  OMEGA 
OMEGA =0MEGA-AL1*XAXIS 
T  =  T 1 

ST  =  SA+SR1 (  I  I 
IF( ST.GE.SY1T  =  0. 

Y1  =  XAXIS 
X1=-L1-YB+(I-1)*AL1 
X2=-L1-YB+I*ALI 
A  =  T*AL  1 

IWY  =  IWY  +  ((0Ml+0MEGA)*Yl*A/2.) 

IWX  =  IWX  +  ( ( 0M1*X1  +  0MEGA*X2» *A/3.  )  +  (  0M1+X2  +  OME  GA*X1 ) *A/6. )  ) 

533  CONTINUE 

C  CALCULATE  OMEGA  SECTION  2 
00  550  1=1, K 
OM 1  =  OMEGA 
0MEGA=0MEGA+AL2*YB 
T  =  T2 
H  =02-1+1 
ST  =  SA  +SR2 ( H ) 

I F { ST.GE.SY) T  =  0. 

Y1  =  XAX  I S  — (  I-  1  >  *AL2 
Y2=XAXIS-I*AL2 
X  1  =  -YB 
A  =  T*AL2 

IWX  =  IWX  +  (0M1+0MEGA )*Xl*A/2. 

IWY  =  IWY  ♦ ( 0M1* Y 1  +0MEGA*Y2 ) *A/3.  +<0M1*Y2  +OMEG A*Y 1 ) *A/6 . 

550  CONTINUE 

C  CALCULATE  LEG  ADDITION 
L  =02- K 
OML  EG=OMEGA 
IFIL.EQ. 0.0)60  TO  1550 
00  1550  1  =  1  , L 
0M1=0MLEG 

0MLEG=0MLEG+AL2*YB 

T=T2 

H=02-K- I +1 
ST=  SA  +  SR2 ( H ) 

IF(ST.GE.SY)T=0.0 
Y 1  =  X  AX  I S-GU AGE- (  I-1)*AL2 
Y2=XAX I S-GUAGF- I  *AL2 
X 1=-YR 
A=T*AL2 

IWX=IWX+(0Ml+0'MLEG)*Xl*A/2. 

I WY= I W Y+ ( OM 1  * Y 1  +0MLEG*Y2 )*A/3.  + I 0M1 * Y2  +  0ML EG* Y 1 ) * A/6. 

1550  CONTINUE 

C  CALCULATE  OMEGA  SECTION  3 
00  560  1=1,03 
OM 1  =  OMEGA 

OMEGA  =  OMEGA  +  ( AL3 * ( GUAG E- X AX  I S )  ) 

T  =  T3 

ST  =  SA  +SR  3  (  I  ) 

IF(ST.GE.SY)T=  .0 
X  1= - Y  B+ (  1-1  )*AL3 
X2  =  - YB+ I *A  L  3 
Y 1  =  ( XAX IS-GUAGE  ) 

A  =  T  * AL3 

IWY  =  IWY  +  ( ( 0Ml+0MEGA)*Yl*A/2. ) 

IWX  =  IWX  +  ( ( OM  1*  X 1  +0MEGA*X2 I *A/3. )  ♦  (  0M1+X2  +  OME GA*X1 ) * A/6 . M 
560  CONTINUE 

CALCULATE  OMEGA  SECTION  4 
CALCULATE  LEG  ADDITION 
OMLEG=OMEGA 
IF(L.EQ.O.O) GO  TO  1560 
001560  1=1, L 
0M1=0MLEG 

OML  EG=0MLFG-AL2  *YB 
T=T2 

H=D2-K- I +1 
ST=  S A+  SR  2 (  H  ) 

IF(ST.GE.SY)T=0.0 
Y1=XAXI S-GU AGE- ( 1-1 ) *AL2 
Y2=XAX IS-GUAGE- I *AL2 
X  1  =  YB 
A=T*AL2 

IWX=IWX+( 0Ml+0MLEG)*Xl*A/2. 

I WY= I WY+ ( QM1 *Y 1 +OMLE  G*Y2 )  *A/3.  ♦ ( 0M1 * Y2 +OML EG*Y 1 ) * A /6. 

1560  CONTINUE 


. 

' 
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AL4=AL2 
DO  5  70  1  =  1 1 K 
0M1=0MEG A 

0MEGA=0MEGA*AL2*YB 

T  =  T2 

H=D2-K*I 

ST  =  SA  +  SR  2 ( H ) 

IF(ST.GE.SY)T=0.0 

Y  1  =  XA  X I S-GUAGE  +(I-1)*AL4 
Y2  =  XAXIS  -  GUAGE  «•  I*AL4 
X  1=YB 

A  =  T  *AL4 

1WY  =  IWY  ♦ ( 0M1 * Y 1  +OMEGA*Y2)*A/3.  +(0M1*Y2  *0MEGA*Y1 )  *A/6. 

IWX  =  IWX  ♦  (0Ml*0MFGA)*Xl*A/2. 

570  CONTINUE 
C  CORNER  SECTION  5 
DO  580  1=1,01 
0M1  =  OMEGA 
0MEGA=0MEGA-AL1*XAXI S 
T  =  T 1 
H=D2-I+1 
ST  =  SA  ♦SRI (HI 
IF(ST.GE.SY)T  =  0. 

Y1=XAX I S 

X  1=YB*(  1-1 )*AL1 
X2  =  Y  B*I  *AL  1 

A  =  T*  AL  1 

IWY  =  IWY  +  ( ( 0M1+0MEGA )*Yl*A/2.  ) 

IWX  =  IWX  +  ( ( OM 1  * X 1  ♦0MEGA*X2) *A/3.  )  +  (  0M1*X2  *■  OME GA* X 1 ) * A/6 .  )  ) 

580  CONTINUE 

WRITE (OUTPUT, 532  IOMEGA, IWY,  IWX 

53  2  FORMAT ( 1  OX, • OMEGA  = ' , FI  0 .4, 8X, • I WY= • , F 8 . 4 , 10X , • I WX= • , F8 . 4 / ) 

534  FORMAT( 10X, • ###########  END  SECTION  ###########  •  /I 

153  2  FORM  AT ( 1  OX, • OML  EG= ' , F 10. 4, 8X , • I W Y= • , F 8 . 4 , 10X , • I WX= • ,F8.4/I 
1534  FORMAT (10X, ' #########END  LEG  ADDITION  ######•/) 

C  CALCULATE  LOCATION  OF  SHEAR  C ENTRE I  TOR S I  ON AL  CONSIDERATIONS) 

XO  =  IWY/IXX 
YO  =-  I WX / I YY 
WRITE(OUTPUT ,581 1X0, YO 

58 1  FORMAT! 1  OX, • X0= • , F 1 0 . 4 , 1  OX , • Y0= • , F 10 . 4/ ) 

C  CALCULATE  WO  SECTION  1  ########### ################## # # ### ##### ###### 

NUM  =  1 

OMEGA  =  0.0 
SN  =  1.0 
YYBAR  =  YBAR-YO 
WRITEIOUTPUT ,600)YYBAR 
600  FORM  AT ( 10X, • YYBAR= •  ,  F12.4/ I 
IF(Y0.GE.YBAR)SN=-1.0 
WO ( 1 1=0.0 
DO  602  I  = 1 , D1 
0MEGA=0 MEGA- AL1 * YYBAR* SN 
NUM  =  NUM  4-  1 
WO(NUM)  =  OMFGA 
ST=S A*SR  1  (  I  I 
IF ( ST .GE.SY) ST=SY 
STP  ESS=ST*F 
Yl= YYBAR 

XI  =  -L1-YB+(2*I-1  )*ALl/2. 

ADIST=X1**2*Y1**2 

BARK  =  BAP,K*2.*STRESS*ADI  ST*AL1*T1 
602  CONTINUE 
C  SFC.T  ION  2 

DO  610  I  =1 , K 
OMEG A=OMEGA  + AL2  * YB 
NUM  =  NUM  +  1 
WO ( NUM I  =  OMEGA 
ST=SA+SR2(  I  I 
IF(ST.GE.SY) ST=SY 
STR  E  SS  =  ST  *  E 

Y  l  =Y Y  BAP- ( 2*I-l)*AL2/2. 

X  1=  YB 

ADIST=X l**2+Yl  * *2 

BARK=8APK+2.*STRFSS*ADIST*AL2*T2 
610  CONTINUE 

WOLL  =OMEGA 

C  CALCULATE  WO  LEG  ADDITION 
OML  EG=OMEGA 
IF(L.FQ.O.O) GO  TO  1610 
DO  1 6 1 0  1=1, L 
0MLEG=0MLEG+AL2*YB 
NUM=NIJM+1 
WO ( NUM) =OMLEG 
ST=SA*SR2(K*I ) 

IF( ST.GF.SY )ST=SY 
STP  E  SS=  S T* F 

Y 1=Y Y  BAR-GU A  GE- (2*1-1 )*AL2/2. 

X  1  =  YB 

ADIST=Y1**2+X1**2 

BARK.  =  BARK+2  .  *STRESS*ADI  ST*AL2*T2 

1610  CONTINUE 
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C  SECTION  3 
SN  =  1.0 

I F ( YYBAR.GE.GUAGEISN  =-1.0 
KK=D3-1 

00  620  I  =  1  , KK 

OMEGA  =  OMEGA  + ( AL 3*ABS ( GUAGE-Y YBAR ) ) *SN 

NUM  =  NUM  +  1 

WO ( NUM I  =  OMEGA 

ST=SA+SR3( I  I 

IE  I ST.GE.SY) ST=SY 

STRESS=ST*E 

Y1=YY  BAR-GUAGE 

X l=YB-( 2*1-1  )  *  AL  3/2 . 

ADIST=Y1**2*X1**2 
BARK=BARK+STRESS*ADI ST*AL3*T3 
620  CONTINUE 

OMEG A=0 MEG A  + AL3 * A BS<  GUAGE-Y YBAR I *SN 
NUM  =  NUM  +  1  +  L 
WO ( NUM )  =  OMEGA 
C  CALCULATE  LEG  ADDITION 
OMLEG=OMEGA 

IF(L.EQ.O.O) GO  TO  1640 
D01 640 1 =1 t  L 
0MLEG=0MLEG-AL2*YB 
NUM=NUM-1 
WO ( NUM) =OMLEG 
1640  CONTINUE 
C  SECTION  4 

NUM=NUM+L 
DO  630  I  =  l,K 
OMEGA=OMEGA  + AL2  *YB 
NUM  =  NUM  +  1 
WO ( NUM )  =  OMEGA 
630  CONTINUE 
C  SECTION  5 

SN  =  1.0 

IFIYO.GE.YBAP )SN=-1.0 
DO  640  1=1, D1 
0MEGA=0MEGA-AL1*YYBAR*SN 
NUM  =  NUM  +  1 
WO ( NUM )  =  OMEGA 
640  CONTINUE 

WR I TE (OUTPUT ,601 ) NUM, OMEGA 
601  FORMAT!  'O', 10X, 'W0(  •  ,13,'  )=•  ,F1 0.4) 

C  CALCULATE  OMEGA  (N) 

C  SECTION  1 
J=1 

OMEGN  =  ,0.0 

WRITE (OUTPUT, 701 ) J, OMEGN 
WN ( 1 ) =0 . 0 
DO  702  I  =  1,01 
T  =  T 1 

ST  =  SA  ♦  SR1III 
IE(ST.GE.SY)T=0.0 

OMEGN=OMEGN+ (W0( J)+WO( J  +  l ) ) *T*AL1 
J  =  J  +  1 
WN(J)  =  OMEGN 
702  CONTINUE 
C  SECTION  2 

DO  710  1=1, D2 
T  =  T2 
H  =  D2-I+1 
ST  =  SA  +  SR2 ( H ) 

I F ( ST.GE.SY) T  =  0.0 

OMEGN= OMEGN* ( WO ( J)+W0(J*1))*T*AL2 
J  =  J  ♦  1 
WN ( J )  =  OMEGN 
710  CONTINUE 
C  SECTION  3 
T  =  T3 

ST=S A+SR3 ( 1 ) 

IF(ST.GE.SY)T=0.0 

OMEGN=OMEGN*(WOLL+WO( J*1 )>*T*AL3 

J  =  J  +  1 
WN  ( J )  =  OMEGN 
DO  720  I  =  2,KK 
T  =  T3 

ST  =  SA  ♦  SR3 !  I ) 

I F ( ST  .GE.SY)T  =  0.0 

OMEGN= OMEGN* (WO ( J )+W0( J+l ) )*T*AL3 
J  =  J  +  1 
WN ( J )  =  OMEGN 
720  CONTINUE 
T  =  T  3 

ST=S A+SR3 ( 03 ) 

IF( ST.GE.SY)T=0.0 

OMEGN= OMEGN* ( W0( J ) +W0( J*L*1))*T*(AL3) 

J= J+L+l 

WN ( J )  =  OMEGN 


C  SECTION  4 

IF(L.EQ.O.O) GO  TO  725 

DO  725  |=1, L 

T=T2 

H=L-I+1 

ST=SA*SR2(  H) 

IF(ST.GE.SY)T=0.0 

OMEGN=OMEGN  *(WO( J ) +WO( J-l ) )*T*AL4 
J  =  J-1 

WN( J)=OMEGN 
725  CONTINUE 
J  =  J  +  L 

DO  730  1=1, K 
T  =  T2 
H=L  ♦  I 

ST  =  SA  +  SR2IH) 

IF ( ST  .GE.SY) T=0.0 
OMEGN=OMEGN-MWO<  J)+WO(  J+l)  )*T*AL4 
J  =  J  +  1 
WN ( J )  =  OMEGN 
730  CONTINUE 
C  SECTION  5 

DO  740  1=1,01 

T  =  T1 

H=01-I*l 

ST  =  SA  +  SR  1  (  HI 

IF(ST.GE.SY)T=0.0 

OMEGN=OMEGN+ (WO ( J J +W0( J +1) )*T*AL 1 

J  =  J  +  1 

WN ( J I  =  OMEGN 

740  CONTINUE 

WR I TE (OUTPUT, 70 1) J, OMEGN 
701  FORMAT < 10X, • WN ( • , 13,  '  )=•  ,F 10.4//) 

AREA=VAREA 

0 ME GAN= OMEGN/ ( AREA *2 . ) 

C  AREA  CALCULATED  ON  REMAINING  ELASTIC  COPE 
WRITE (OUTPUT , 74 1 ) OMEGA N, ARE A 

741  FORM  AT ( 10X , • OMEGAN=* , F 10. 4 , • ARE A= ' ,F10.4//I 
C  CALCULATE  WNN 

C  WNN=W( N) /AREA 
SUM=0.0 
DO  802  1=1, J 
COR=WO ( I ) 

AWNN=OMEGAN-COR 
WNN ( I ) = AWNN 
SUM=  SUM  +  WNN ( I ) 

802  CONTINUE 

WRITE (OUTPUT, 801 ) I , WNN ( I ) 

801  FORMAT!  10XK'WNN(  •  ,  I  3  ,  •  )  =  •  ,  F  10 . 4/ ) 

C  IF  DISTRIBUTION  OF  PREVIOUS  PRORERT I ES  OK  SUM=0.0 
WR I TE (OUTPUT ,804 ) SUM 
804  FORMAT ( 10X, • SUM=* ,F 10. 6/ 1 
C  CALCULATE  WARPING  MOMENT  OF  INERTIA(IW) 

C  SECTION  1 

AA  I  W  =  0 . 0 
J=1 

DO  830  1=1,01 
T=T  1 

ST=  S  A+SR 1 ( I > 

IF(ST.GE.SY)T=0.0 

AAIW=AAIW-M  (WNN(  J) **2 )  ♦  (  WNN  (  J  )  *WNN  (  J  + 1  >  )  ♦  (  WNN  (  J  ♦  1  )  **2  ))  *  T*  AL  1 
J  =  J+1 

AIW(J)=AAIW 

830  CONTINUE 
C  SECTION  2 
T=T2 

D0840  1=1, 02 
T  =  T2 

H=  02- I  +  1 
ST=SA-*-SR2(  H  I 
I F ( ST.GE . SY 1 T  =  0.0 

AAIW=AAIW-M  ( WNN ( J ) **2 )  ♦  ( WNN ( J ) *WNN ( J+ 1 1 >♦ ( WNN < J ♦ 1 ) **2 ) ) * T* AL 2 
J  =  J  +  1 

A  I W ( J)=AAIW 
840  CONTINUE 
C  SECTION  3 
T  =  T  3 

ST=SA*SR3( 1 ) 

IF(ST.GE.SY)T=0.0 

M=J-L  .  .  , 

AA  I  W=  AA  I  W  +  (  ( WNN(  M  )  **2  )  *■  (  WNN  (  M  )  *WNN  (  J  ♦  1 )  )  *  ( WNN (  J  +  1  )  **2  )  )  *T* AL  3 

J  =  J+1 

AIW( J»=AAIW 

DO  850  1=2, KK 
T=T3 

ST=S A+SR3 ( I ) 

IF(ST.GE.SY1T=0.0  . 

A  A  I W  =  AA I W+(  ( WNN ( J ) **2 )  +  ( WNN ( J ) *WNN ( J ♦ 1 1 ) ♦ ( WNN ( J ♦ 1  )  **2 I  I  * T * AL 3 

J  =  J  +  1 

AIW(J)=AAIW 

850  CONTINUE 
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T=T  3 

ST=SA+SR3(03 ) 

IF! ST.GE.SY)T=0.0 
N= J+L+l 

AAIW=AAI W+( (WNN! J)**2 I  ♦  ( WNN (J ) *WNN ( N ) ) ♦ < WNN ( N ) **2 I) *T *AL 3 
J  =  J  +  1 

AIW! J)=AAIW 
C  SECTION  4 

DO  860  1=1, D2 
T=T2 

H=D2-K+1 
ST=SA*SR2( H) 

IE(ST.GE.SY)T=0.0 

AAIW=AAIW-H  (WNN( J)**2  )  +  IWNNIJ )  *WNN  (  J  1 1 )  ♦ ( WNN ( J+ 1  1**2  )  )  *T*AL  2 
J  =  J+1 

AIW( J I =A A I  W 
860  CONTINUE 
C  SECTION  5 

DO  870  1=1, D1 
T=T  1 

H=D 1- I *1 
ST  =  S  A+SR 1(H) 

IF( ST.GE,SY)T=0.0 

AAIW=AAIW*( (WNN! J)**2)  ♦  (  WNN  ( J  )  *WNN  (  J  ♦  I )  )  ♦  ( WNN  (  J  ♦  1  )  **2  )  )  *T  *AL  l 
J=  J  +  l 

AIW( J)=AAIW 

870  CONTINUE 
WRITEIOUTPUT ,831  )J,AAIW 

831  FORMAT! 'O' ,10X, • IW( • ,  13, • )=' , FI  0.4) 

A  I WW=A A  I W/3  . 

WRITEIOUTPUT ,871  IAIWW 

871  FORMAT! 1  OX , •  IW  OF  SECT  1 0N= • , FI  0 . 4/ / ) 

C  CALCULATE  VALUES  FOR  8UCKLING  EQUATIONS 
C  NOTE  (R0**2 I=BARK/LOAD 

XA=E*IYY*3. 14159*3. 14159 
XB=LOAO 

XC=  E*  A  I WW*  3 . 14159*3. 14159 

XD=G*AKT 

XE=RARK 

XF=(L0A0*Y0)**2 
WRI TE l OUTPUT ,902 ) XA , XB , XC 
WRITE (OUTPUT, 902) XD.XE.XF 
902  FORMAT! 10X,F14.4,3X, F14.4,3X,F14.4/I 
C  SET  UP  QUADRATIC  FOR  SOLUTION  IN  TERMS  OF  A  CRITICAL  LENGTH 
C  QUADRATIC  (  AA(L**4)  ♦  RB(L**2)  +  CC  =  0.0 
AA=-XB*XD-XB*XE-XF 
BB=XA*XD*XA*XE-XB*XC 
CC=XA*XC 

WRITE (OUT PUT, 904 ) AA,BB,CC 

904  FORM  AT! 10X, ' A=*  , F14.4,3X, *B= • ,F14. 4, 3X , *C=' ,F14.4/ ) 
XAC=(BB*BB-4.*AA*CC) **0.5 
XX1=(-BB+XAC)/(2.*AA» 

XX2=(-BB-XAC)/(2.*AA) 

WRITE(OUTPUT ,905)XX1,XX2 

905  FORM  AT (10X,F14.4,8X,F14.4//) 

C  CHECK  FOR  AN  IMAGINARY  SOLUTION 
C  XXI  AND  XX2  REPRESENT  SOLUTIONS 

IF! XX1.GE.0.0)LL1=XX1**0. 5 
IF(XX2.GE.0.0)LL2=XX2**0.5 
WR I TE! OUTPUT, 900  ILL 1,LL2 

900  FORMAT! 1 OX , • CR I T I CA L  L ENGTH= • , F 1 0 . 4 , 8X , «0R= ' , F 10 . 4 ,  '  I NC HE S • / I 
LOAD=0 . 0 
IXX=0.0 
I YY=0. 0 
A  =  0 . 0 
AD=0 . 0 
GO  TO  10 
999  STOP 

C  SAMPLE  DATA  LAYOUT 

C  MODULUS  OF  ELASTICITY,  YIELD  STRESS,  SHEAR  MODULUS 

C  SECTION  GEOMETRY  AND  NUMBER  OF  DIVISIONS  FOR  EACH  LEG, GAUGE  LOCATION 
C  RESIDUAL  STRAIN  FOR  EACH  DIVISION 

C  APPLIED  STRAIN  LEVEL ( I NCREMENT  UNTIL  SECTION  FULLY  YIELDEDI 
END 

ADCONS  FOR  EXTERNAL  REFERENCES 
******  eno  nF  COMPILATION  ****** 


